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Abstract: We propose a model for a power-counting renormalizable field theory 
living in a fractal spacetime. The action is Lorentz covariant and equipped with a 
Stieltjes measure. The system flows, even in a classical sense, from an ultraviolet 
regime where spacetime has Hausdorff dimension 2 to an infrared limit coinciding 
with a standard D-dimensional field theory. We discuss the properties of a scalar field 
model at classical and quantum level. Classically, the field lives on a fractal which 
exchanges energy-momentum with the bulk of integer topological dimension D. Al- 
though an observer experiences dissipation, the total energy-momentum is conserved. 
The field spectrum is a continuum of massive modes. The gravitational sector and 
Einstein equations are discussed in detail, also on cosmological backgrounds. We find 
ultraviolet cosmological solutions and comment on their implications for the early 
universe. 
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1. Introduction 

The search for a consistent theory of quantum gravity is one of the main issues in the 
present agenda of theoretical physics. Beside major efforts such as string theory and 
loop quantum gravity, other independent lines of investigation have recently attracted 
some attention. Among these, Hofava-Lifshitz (HL) gravity [jl], Q is a proposal 
for a power-counting renormalizable model 0, Q which is not Lorentz invariant. 
Coordinates scale anisotropically, i.e., [t] = —z and [x*] = — 1 in momentum units, 
where ^ > 3 is a critical exponent typically fixed at 2; = 3. Because of this, the total 
action can be engineered so that the effective Newton constant becomes dimensionless 
in the ultraviolet (UV) and higher-order spatial derivatives improve the short-scale 
behaviour of particle propagators. Due to the presence of relevant operators, the 
system is conjectured to flow from the UV fixed point to an infrared (IR) fixed point 
where, effectively, Lorentz and diffeomorphism invariance is restored at classical level. 

Another property of HL gravity stemming from the running of the couplings 
effective dimension is that the spectral dimension (is [@, |, ^ H at short scales is 
(is ~ 2 0. This is in intriguing accordance with other proposals for quantum gravity 
such as causal dynamical triangulations , asymptotically safe gravity |Ty] and spin- 



foam models |jTT|] (see also [T^). Systems whose effective dimensionality changes with 
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the scale can show fractal behaviour, even if they are defined on a smooth manifold.^ 
All the above examples incarnate the popular notion that "the Universe is fractal" 
at quantum scales. 

Despite the beautiful physics emerging from the HL picture, inspired by critical 
and condensed-matter systems, it potentially suffers from at least one major problem. 
Lorentz invariance, one of the best constrained symmetries of Nature, is surrendered 
at fundamental level. As argued on general grounds |T^, ITBl, even if deviations 



from Lorentz invariance are classically negligible, loop corrections to the propagator 
of fields lead to violations several orders of magnitude larger than the tree-level 
estimate, unless the bare parameters of the model are fine tuned. This expectation 



1^ is indeed fulfilled for Lifshitz-type scalar models fl^. Although supersymmetry 



might relax the fine tuning [0, the present version of HL gravity is clearly under 
strong pressure, also for other independent reasons. 

Motivated by the virtues and problems of HL gravity, it is the purpose of this 
paper to formulate an effective quantum field theory with two key features. The 
first is that power-counting renormalizability is obtained when the fractal behaviour 
is realized at structural level, i.e., when it is implemented in the very definition of 
the action rather than as an effective property. In other words, we will require not 



only the spectral dimension of spacetime, but also its UV Hausdorff dimension [19 



(which will coincide with ds in our case) to be cin ~ 2. Secondly, we wish to maintain 
Lorentz invariance. 

Therefore, this proposal is (a) defined on a fractal (in a sense made precise 
below), (b) Lorentz invariant, (c) power-counting renormalizable, (d) UV finite with 
no ghost or other obvious instabilities, and (e) causal. A condensed overview of the 
model was given in |^ . 



Some of the ingredients we shall use are similar to those found in other recipes 
(e.g., scalar-tensor theories or models with fractional operators). Their present mix- 
ing, however, will hopefully give fresh insight into some aspects of quantum gravity. 
For example, a running cosmological constant naturally emerges from geometry as a 
consequence of a deformation of the Poincare algebra. 

The plan of the paper is the following. The main idea is introduced in section 
1^. With particular reference to a scalar field theory, a dimensional analysis of the 
coupling constants is given in section |2.1| . Section ^ is devoted to a scalar field on a 
Minkowski fractal: its classical equation of motion and dynamics are presented in sec- 
tion 3A, where the Hamiltonian formalism is shown to admit both a dissipative and 
conservative interpretation. The causal propagator of the free field in configuration 
space is calculated in section p.2| , while its Fourier-Stieltjes transform in momentum 
space is discussed in section p?3| . We outline the gravitational sector in section |[ 

^HL gravity with detailed balance possesses a natural fractal structure also because of the ap- 
pearance of fractional pseudo-differential operators This version of the theory, however, 
seems to be unviable. 



-2- 



Einstein and cosmological equations are derived in sections iA and i2, respectively, 
where cosmological solutions are found and analyzed. Section ^ contains concluding 
remarks and a discussion on open issues and future developments. 



2. Fractal universe 

In HL gravity, one requires that time and space coordinates scale anisotropically. On 
one hand, this leads to a running scaling dimension of the couplings and an effective 
two-dimensional phase in the UV. On the other hand, anisotropic scaling gives rise 
to higher-order spatial operators and a non-Lorentz-invariant action. It turns out 
that we can achieve the first result (and avoid the second) by maintaining isotropic 
scaling, 

[x^] = -l, /i = 0,l,. (2.1) 
while replacing the standard measure with a nontrivial Stieltjes measure, 

d^x dg{x) , [g] = -Da ^ -D . (2.2) 

Here D is the topological (positive integer) dimension of embedding (abstract) space- 
time and a > is a parameter.^ What kind of measure can we choose? A two- 
dimensional small-scale structure is a desirable feature of renormalizable spacetime 
models of quantum gravity, and the most naive way to obtain it is to let the effective 
dimensionality of the universe to change at different scales. A simple realization of 
this feature is via fractional calculus and the definition of a fractional action.^ 



To begin with, we quote the following results in classical mechanics. In 
empirical evidence was given that the Hausdorff dimension of a random process 
(Brownian motion) described by a fractional differintegral is proportional to the 
order a of the differintegral; the same relation holds for deterministic fractals, and in 
general the fractional differintegration of a curve changes its Hausdorff dimension as 



dn ^ dn + a (see also ||2^). Moreover, integrals on net fractals can be approximated 
by the left-sided Riemann-Liouville fractional integral of a function L{t) [EB, 12^ 



^The Hausdorff dimension of a set is greater than or equal to its topological dimension but 
the situation one has in mind here is a physical spacetime (the fractal) embedded in an ambient 
D-dimensional manifold M . All physics takes place in the fractal and there are no observers in the 
"bulk" Ai. Given this picture, one can interpret the present model as "diffusion of spacetime" in 
an embedding manifold. 

^Another route, which we shall not follow here, is to define particle physics directly on a fractal 
set with general Borel probability measure g. This was done in (and [^2[ ^ on Sierpinski 
carpets) for a quantum field theory on sets with Hausdorff dimension 4 — e very close to 4. The 
model in |^ has many aspects of dimensional regularization ^ , one difference being that the 
parameter e is taken to be physical. 
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31, 32 , 



dg{t) L{t) oc QliL{t) 



g{t) 



-i- f dt{i-tr"'L{t) 
r(a) io 

r(a + l) ' 



(2.3) 
(2.4) 
(2.5) 



where t is fixed and the order a is (related to) the Hausdorff dimension of the set 

^ is vahd for large Laplace momenta and can 



2^, |3^ . The approximation in eq. 
be refined to better describe the full structure of the Borel measure g characterizing 
the fractal set. In the latter case, integration on the set is approximated by a sum 
of fractional integrals ppj] . 

Different values of < a < 1 mediate between full-memory (a = 1) and Markov 
processes (a = 0), and in fact a roughly corresponds to the fraction of states pre- 
served at a given time t during the evolution of the system [|^, ^ Applications 
of fractional integrals range from statistics, diffusing or dissipative processes with 
residual memory [^, such as weather and stochastic financial models [0, to sys- 



tem modeling and control in engineering [p5 |. 
Noticing that a change of variables t i 



t transforms eq. (|2.4| ) into the form 



r(a) 



dte-'L{t-t) 



(2.6) 



the Riemann-Liouville integral can be mapped onto a Weyl integral in the limit 
i — )■ +00. The limit is formal if the Lagrangian L in eq. ( p.6| ) is not autonomous. 
We assume otherwise, so that limt_^oo -^(^~ i) = L[q{t),q{t)]. 

This form will be the most convenient for defining a Stieltjes field theory action. 
When i — )■ +oo, eq. ( |2.6| ) is proportional to the usual formula in a dimensions em- 
ployed in dimensional regularization. After constructing a "fractional phase space" 
1^, ^ , this analogy confirms the interpretation of the order of the fractional 



integral as the Hausdorff dimension of the underlying fractal ^8 



All the above results in one dimension can be easily generalized to a D-dimensional 
Euclidean space (e.g., |^ 0), thus opening a possibility of applications in space- 
time. We entertain the possibility of formulating a scalar field theory with Stieltjes 
action, for the purpose of controlling its properties in the ultraviolet.^ The scalar 
field model is interesting in its own right but also as a simple example whereon to 



^Introductions on the Lebesgue-Stieltjes integral can be found in ji^, |4|, A neat geometrical 
interpretation of the Riemann-Stieltjes one-dimensional integral as the projected "shadow of a 
fence" is given in 47]. Projection is a tool sometimes employed to determine the Hausdorff 
dimension of a fractal llEl. 
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work out the physics. After that, we shall explore the gravitational sector. In D 
dimensions, we consider the action 

S = [ dg{x)C{<l),d^<P), (2.7) 
Jm 

where C is the Lagrangian density of the scalar field (f){x) and 

D-l 

dg{x) = l[f(^){x)dx^ (2.8) 

At=0 

is some multi-dimensional Lebesgue-Stieltjes measure (actually a Lebesgue measure, 
if g is absolutely continuous, which we assume to be the case) generalizing the triv- 
ial Z)- dimensional measure d^x. We denote with {Ai,g) the metric spacetime Ai 
equipped with measure g. We shall consider the situation where is a manifold 
but this may not be the case in general. 

Equation ( p^.7| ) resembles a field theory with a dilaton or conformal rescaling 
V = Yl^ f{pL) of the Minkowski determinant. As in these other models, one will 
obtain an extra friction term in the equation of motion, although the physics will 
be radically different both at microscopic and macroscopic level. This is because 
the measure weight must scale in a certain way, while dilaton solutions in effective 
actions of string theory typically enjoy much more freedom. We should stress at least 
two more reasons why the present model is not just an exotic reformulation of dilaton 
scenarios.^ First, the dilaton of string theory couples differently in different sectors, 
thus leading to a violation of the strong equivalence principle; in our case, the scalar 
field V is still of geometric origin but appears as a global rescaling. Second, a change 
in the measure is accompanied by a new definition of functional variations and Dirac 
distributions, in turn leading to an unfamiliar propagator and the deformation of the 
Poincare algebra. 

If g is not invariant under the Lorentz group S0{D — 1^ 1), the equipped manifold 
is not isotropic even if M. is the Minkowski flat manifold. If only global Poincare 
invariance is broken, as it typically happens in fractals, (A^, g) is not homogeneous. 
Since we wish the Lorentz group SO{D — 1,1) to be part of the symmetry group of 
the action, the Lagrangian density £ and the D weights /(^) must be Lorentz scalars 
separately. The former can be taken to be the usual scalar field Lagrangian, 

C = ~^d,m-V{<P), (2.9) 

where is a potential and contraction of Lorentz indices is done via the Minkowski 
metric r]^i, = ( — h ■ ■ ■ +)fj.u- As for the Stieltjes measure, we make the spacetime 
isotropic choice 

/m = /, /i = 0,l,...,Z}-l. (2.10) 

^ These reasons do not forbid the fractal model to admit also a "dilaton-like" reformulation (see 
below). 
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This should eventually correspond to a fractal in time and space. There are many 
other Ansdtze, for instance an isotropic nontrivial measure 



/{o) = l, /{.) = /, i = l,...,D-l, (2.11) 
or an anisotropic measure of the form 

These measures will correspond to different dynamics but, by construction, to the 
same UV Hausdorff dimension c/h ~ 2. We take eq. ( p.lO| ) with 

v^f"", H = D(l-a), (2.13) 

which generalizes eq. ( |2.6| ).^ The scalar field action reads 

'1 



S = - I d^xv 



9^09^0 + 1^(0) 



(2.14) 



We now pause and discuss the interpretation of the measure. Classically, one can 
boost solutions of the equation of motion to a Lorentz frame where v = f (x) (space- 
like fractal) or f = v{t) (timelike fractal). These two cases will lead to different 
classical physics but at quantum level all configurations should be taken into ac- 
count, so there is no quantum analogue of space- or timelike fractals.^ In any case, 
we shall see that the theory on the Da -dimensional fractal is expected to be dissipa- 
tive, i.e., nonunitary. This conclusion is in line with the known results of fractional 
mechanical systems. Fortunately this will not be a problem because, from the point 
of view of the manifold Ai with D topological dimensions, energy and momentum 
are indeed conserved. 



2.1 Renormalization 

The scaling dimension of is 



Da -2 , , 

(2.15) 



6 



V — v{x) is a coordinate-dependent Lorentz scalar. An alternative generalization of non- 
relativistic fractals might have been to choose v to be also metric dependent, e.g., v = 
igfii/X^X'^l ^ for some vector x^. In this case, however, one does not obtain a consistent set of 
equations of motion. 

^We have seen that, approximately, fractional integrals can represent systems living on a certain 
class of fractal sets. Since we will assume the existence of a nontrivial renormalization group flow 
entailing integrals of different orders ai,a2, ■ ■ ■ , the complete all-scale picture beyond classical level 
will not be a fractal with scale-independent Hausdorff dimension but a multifractal. For this and 
the reason stated in the text, the fractal interpretation of fractional integrals is more involved in 
the quantum theory. To our purposes it is not necessary to stick with it, although we shall do so 
with a slight abuse of terminology. We shall see that the spectral dimension of the universe changes 
in a precise, a-dependent way, thus justifying the term "fractal" a posteriori. 
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which is zero if, and only if, 

a = ^. (2.16) 

Then a = 1/2 in four dimensions. This value can change for the other measures 
defined in eqs. ( |2.10| ) and ( ^.12 ). In the example ( ^.12| ) with [/] = 1 — a, one has 



a = 2/{D — i), and in order for the integral to be properly fractional (rather than 
multiple) it must he i < D — 3. In four dimensions, there can be at most one ordinary 
direction, li i = 1, then a = 2/3. i = corresponds to eq. (|2.13|) , which we shall 
adopt from now on. 

Let the scalar field potential be polynomial. 



N 



v = j2^nr, (2.17) 



n=0 

and let N be the highest (positive) power. The coupling has engineering dimen- 
sion 

, , N(Da-2) 

M = Da-^^ ^. (2.18) 

For the theory to be power-counting renormalizable [cttv] > 0, implying 

N<-^ if a>^, (2.19) 
Da — 2 U 

N<+oo if a<^. (2.20) 

When a = 1, one gets the standard results [0] = {D — 2)/2, N < 2D/{D — 2); in four 
dimensions, the 0^ theory is renormalizable. In two dimensions, is unconstrained. 

These considerations lead us to try to have the parameters run from an ultraviolet 
nontrivial fixed point where a = 2/Z? to an infrared fixed point where, effectively, 
d = aiR. The dimension of spacetime is well constrained to be 4 from particle physics 
to cosmological scales and starting at least from the last scattering era ^ ^ . 

Therefore, = 1 if D = 4. To actually realize this particular flow, one should add 
relevant operators to the action corresponding to terms with trivial measure weight. 
The total scalar action then is 



S 



/.-4„£ + M''<.-".£]. (2,21, 



where M is a constant mass term ([M] = 1) and C is C, eq. ( p.9|) , with all different 
bare couplings (ct„ — )■ 5"„). We symbolically represent this modification of the action 
as 

v{x) -> v{x) + M^(^-°) . (2.22) 

The constant term is anyway expected in the most general Lorentz-invariant defini- 
tion of the measure weight. 
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A fractal structure must shortly evolve to a smooth configuration. Oscillations 
of neutral B mesons can constrain the typical UV mass scale to be larger than about 

m 

M > 300 ^ 400 GeV. (2.23) 

Here we do not attempt to place constraints on this scale with other high-energy 
observations. 



As one will see in section |3.3| , convergence of the Feynman diagrams is better 
than in four dimensions, as one can check by looking at the superficial degree of 
divergence, which is the same as for a Da- dimensional theory [^. In the case of 
gravity, in fact, the usual configuration-space results in 2 + e dimensions should apply 
near the UV fixed point |3|, HIIHI, [S^, |58|, |59 . 



Needless to say, the above construction and remarks fall short of demonstrating 
the existence and effectiveness of such a fiow, which should be verified by explicit 
calculations. Our attitude will be to introduce the model and first see its character- 
istic features and possible advantages, leaving the issue of actual renormalizability 
for the future. 

At any rate, classically the system will fiow from a lower- dimensional fractal 
configuration to a smooth D-dimensional one. This is clear from the definition ( p.22| ) 
of the measure weight and its scaling properties when a < 1. At small space-time 
scales, the weight v ~ dominates over the constant term, while at large 

scales it is negligible. This is true simply by construction, and independently from 
renormalization issues. 

Therefore, the phenomenological valence of the model is guaranteed, at least. 
In our framework both the Newton's coupling and the cosmological constant will 
vary with time already at classical level. On one hand, in minisuperspace models 
motivated by other approaches to quantum gravity, the running of the couplings 
can be implemented at the level of the equations of motion, thus obtaining a high- 
energy "improved" dynamics. This strategy is adopted, for instance, in the Planck- 
ian cosmology of asymptotically safe gravity thanks to its renormalization properties 
[pO| , |61| , |6^ ]. On the other hand, a phenomenological time- varying dimension can 



be considered for constraining the transition scale from fractal to four- dimensional 
physics [^. The couplings running is then also obtained in fractal- related cosmo- 
logical toy models with variable dimension [p^.^ 

3. Scalar field theory 

3.1 Equation of motion and Hamiltonian 

The Euler-Lagrange and Hamilton equations of classical mechanical systems with 



*A11 these scenarios differ in philosophy with respect to [Q, where the spacetime dimension is 
promoted to a dynamical field. 
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(absolutely continuous) Stieltjes measure have been discussed in |6^, ^ ^ in the 



one-dimensional case and ||68| , |69| [fO|| in many dimensions. The Euler-Lagrange 



equation of scalar field theory can be found in |]7I|. 

We can easily adapt the same procedure in our case. From now on we consider 
only the UV part of the action, setting M = 0. Any result in the infrared can then 
be obtained by going to the effective limit a — t- 1. 

The metric space is equipped with a nontrivial measure and caution should be 
exercised when performing functional variations. For instance, the correct Dirac 
distribution is 

1 = / dQ{x)5\l'\x) , (3.i; 



as was also noticed in ■ Invariance of the action under the infinitesimal shift 



+ (50 

yields the equation of motion (for a generic weight v 





dC 



daV d 



dC 



From eq. ( |2.14| ) we get 



□0 + -^^'^0 -v = 0, 



(3.2) 



(3.3) 



(3.4) 



where □ = d^d'^ and a prime denotes differentiation with respect to 0. 

The above friction term is characteristic of dissipative systems and one would 
expect energy not to be conserved. In fact, the Hamiltonian is no longer an integral 
of motion. Let us define the momentum 



6S 



(3.5) 



where dots indicate (total) derivatives with respect to time and we have taken 
eq. ( |3.1| ) into account. 

Defining the Lagrangian 

L = I d'^vC , (3.6) 



the Hamiltonian is 



H = J (ix t> (7r00) — L 

= j dxv{^7ll + ldi(l)d'^+V) 



(3.7) 



(3. 
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where = dx^ . . . dx^ ^. The definition of the equal-time Poisson brackets (time 
dependence imphcit) 



{A(x),i?(x')}. = J dgiy) 
yields the Hamilton equations 



M(x) 55(x') M(x) 55(x') 



(3.9) 



{^,H%, (3.10) 



V 



T^c/) = {T^fl}, H}^ - -TT^ , (3.11) 



equivalent to eqs. (|3.5| ) and ( |3.4|) , respectively. Therefore, time evolution of an ob- 
servable 0{(p,7iip,x) is 

d = dtO + {O, H}, - \^ . (3.12) 



Equations ( |3.11| ) and ( ^.121 ) signal dissipation. Nonconservation of the Hamiltonian 



is expected from the definition of the Lagrangian (|3.6| ). The measure factor v is both 
time and space dependent, so it was not possible to factorize it in order to write the 
action as the correct Stieltjes time integral of L. However, one can exploit Lorentz 
invariance and pick a frame where v = v{\x.\). In this frame H would be restored as 
the generator of time translations; for a timelike fractal it is not possible to pick this 
frame and there is always energy dissipation. Conversely, physical momentum will 
be dissipated in a spacelike fractal. 

Is there a problem with that? Whatever the choice of classical fractal model, 
one would have to face the issue of unitarity at quantum level. Moreover, we need 
a physical interpretation of dissipation. It turns out that the latter helps to address 
the above concern. 

Consider a (-D — 1) -dimensional box of size / and spatial volume l^"^. At the 
scale I, particles live effectively in Da spacetime dimensions. If a = 1, they occupy 
the whole phase space in the box. Otherwise, they must dissipate energy, since the 
energy of the configuration filling the entire topological volume is different from that 
of a configuration limited to the effective Da-dimensional world. The total energy of 
the system S in D topological dimensions is conserved, but the energy H measured 
by a Da-dimensional observer is not (eq. (|3.12|) with O = T-Lv and integrated in 
space): 



H = dtH- J dxwTr; 

= - / rfx{;£. (3.13) 
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In fact, a conserved quantity is 



S{t) = H{t) + A{t) 
A(t) = / dt I dy.vC. 



(3.14) 
(3.15) 



which we are going to obtain also from Noether's theorem (see ^ for a com- 
putation in classical mechanics). It should be easy to see that A can be regarded as 



the complementary function of initialized fractional calculus |73|, |7^, 0. Phys- 
ically, it is a running cosmological constant of purely geometric origin. For this 
reason, dissipation might eventually prove to be an asset rather than a liability of 
the theory. 

Unlike standard scalar field theory, the Noether current associated with the 
usual Lagrangian continuous symmetries is not covariantly conserved. On the other 
hand, one can easily find generalized conserved currents. Take a generic infinitesimal 
transformation of the field, eq. (|3.2|) , and coordinates, ^ + 5x^. We consider 
symmetries of the autonomous Lagrangian density C and define "quasi invariance" 



of the action as done in |67, 71, 72|. Then 5C is a total divergence, 



dj^ 
dxf^ 



(3.16) 



Combining this equation with 
6C = ^5<h + 



dC 



:Kd,<P) 



dC 



d dC 



\ d(j) dx^^ dd^ 



+ 



d 



dxi^ V dd, 



dC 



one obtains on shell 



d 



dx^" 



dC 

dd,d 



+ J^d^v = 



Choose a coordinate translation 



x^" = x^' + Sx" , 5x'' 

so that 5(f) = a''d^(j) and SC = a'^dC/dx^. Then 

d 



(3.17) 



(3.18) 



(3.19) 



(3.20) 
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where 

Tt = b^,C + d^<^d,<^ (3.21) 
is the usual energy-momentum tensor. Integrating eq. (|3.20|) in space, one gets 



Py + / dx(9^t;£ = 



where 



(3.22) 



(3.23) 



The 1/ = component yields Pq = H and eq. (|3.13 ) follows suit. The v = i component 
gives the physical momentum 



and its conservation law 



Pi + ki = Pi+ d^divC = 



(3.24) 



(3.25) 



Pi generates spatial translations in the field but not in its conjugate momentum, 
since the covariant counterparts of eqs. ( ^.1U| ) and ( |3.11| ) are 



di,(p = {0, Pf,}v , 
d^j.iT^ = {rc^, Pfj,}^ ^TT^ 



The covariant version of eq. (|3.12 

dO 



IS 



yielding 



dx^^ 



dP, 



dxf^ 



d,0 + {0,P,%-^n^^ 



(3.26) 
(3.27) 

(3.28) 
(3.29) 



The (/i, u) components of this equation offer consistency checks and new commutation 



relations. Component (0,0) corresponds to eq. (|3.13|) , while gives 

{P„Pj}, = 0. (3.30) 



Components {0,i) and (i,0) and eq. ( |3.25| ) are consistent among each other if, and 
only if. 



{H,P,}, 



c/x diV C , 



(3.31) 



finally showing that the Poincare algebra is now noncommutative unless v is only 
time dependent (timelike fractal). One can check that also the Lorentz algebra is 
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deformed. Consider the Noether current associated with boost /rotation transforma- 
tions 



50 = a^^ix'^d'' - , (3.32) 

6C = a^aix^d" - = a^,9^(^^'"x^£ - g^'^x'^C) = d^J^ . (3.33) 



Substituting in eq. (|3.18|) , we get 

d 



{vM>""') + C^x^d^v - x^d^v) = , (3.34) 



dxf^ 
where 

M^""" = x'^T^" - x^T^"" . (3.35) 

The algebra of the Lorentz generators J''"" = J d'x.vM^'^" ^ which includes the Lorentz 
boosts = and the angular momenta = ^e^j^J-'^ jk the Levi-Civita 
symbol), is deformed because of the nontrivial weight v. 

We have ended up with a classical field theory living on a Z^a-dimensional fractal 
breaking Poincare invariance. Depending whether the fractal is "timelike" or "space- 
like," there will be a momentum or energy transfer between the world-fractal and 
the D-dimensional embedding. 

The D-dimensional side of the picture can be actually made more precise. So 
far we have interpreted the function v in the action (|2.14|) as (the derivative of) a 
Stieltjes measure defined on a fractal of Hausdorff dimension Da. Of course one 
can regard it as a "dilaton" field coupled with the Lagrangian density C living on a 
D-dimensional manifold. Then it is natural to consider the usual 5 of Dirac, 

d^a;5(^)(a;). (3.36) 

Consequently, the momentum conjugate to the field cj) is 

vr^ = vn^ . (3.37) 
It is easy to convince oneself that the Poisson brackets are the usual 

(3.38) 

and that all v dependence disappears in the D-momentum 

P^ = P^((/.,7r,). (3.39) 

Poincare (and Lorentz) invariance is preserved. This completes the proof that, at 
least at classical level, dissipation occurs relatively between parts of a conservative 
system. Quantization would follow through, although an UV observer would experi- 
ence an effective probability flow from or into his world-fractal (see also 
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To summarize, because the universe is associated with a "fractal" structure one 
typically expect to have breaking of Lorentz invariance at small scales. The geom- 
etry of the problem is not standard and this modifies the usual definition of Dirac 
distribution and Poisson brackets. As a result, from the point of view of an observer 
living in the fractal and measuring geometry with weight v, Lorentz and translation 
invariance are broken inasmuch as the action itself is Lorentz and Poincare invariant, 
but the algebra of the Poincare group is deformed. In other words, when talking 
about fractals embedded in Minkowski spacetime we mean the geometries defined 
by the deformed Poincare group. However, from the point of view of the ambient 
D-dimensional manifold Poisson brackets and functional variations no longer feature 
the nontrivial measure weight, which is now regarded as an independent matter field 
rescaling. In that case, the full Poincare group is preserved. On the other hand, the 
properties of nonrelativistic fractals are more intuitive: they break translation in- 
variance because the measure weight introduces explicit coordinate dependence and 
the system is not autonomous. 

3.2 Propagator: configuration space 

The theory is Lorentz invariant, ghost free and causal at all scales. We can check 
this explicitly by computing the causal propagator G{x), which is (proportional to) 
the propagator in two and D dimensions at the UV and IR fixed points, respectively. 
As usual, we define the vacuum-to- vacuum amplitude (partition function) 

Z[J] = J[d(j)]e'^'^^^^+'l'^\ (3.40) 

where J is a source and we have already integrated out momenta. Integration by 
parts in the exponent allows us to write the Lagrangian density for a free field as 

jC^^(P(d+ - m^^ = ^(PC(j) . (3.41) 

The propagator is the Green function solving 

CG'(x) = 5f)(x). (3.42) 

By virtue of Lorentz covariance, the Green function G must depend only on the 
Lorentz interval 

= x^x^ = Xix' - t^ , (3.43) 

where t = x^ and i = l,...,D — 1. In particular, v = v{s) with the correct scaling 
property is v{s) ~ This definition guarantees reality of the measure and 

avoids problems with unitarity (in particular, the action is real). 

One might be worried that the measure blows up on the light cone, but this is 
an integrable singularity: the check that d^x v < oo is done on a compact set A 
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(for instance, a D-hall of radius R) and in Euclidean signature. Below we will also 
see that there is nothing pathological in the propagator on the light-cone, even if v 
is singular in s = 0. Measures describing fractals may be very irregular and it would 
be worth investigating the physical interpretation of their singularities, especially in 
Lorentzian signature. For the time being, we notice that one can also define the 
measure weight to be 

CI |D(q-1) /q 

This definition may be in contrast with the original assumption that the measure be 
absolutely continuous. This would mean that the simplified model with an overall 
measure weight does not come from a most general fractal model with Lebesgue- 
Stieltjes measure g: in other words, dg{x) ^ v{x) dPx. This is not a problem for two 
reasons. On one hand, the simple model with measure weight v is still able to capture 
much of the physics of the general Stieltjes model, by virtue of the scaling argument; 
as a matter of fact, one could even take measures which are Lebesgue-Stieltjes only 
asymptotically, and yet obtain a modelization of a fractal quantum field theory in 
certain regimes. On the other hand, one can devise other measure profiles with the 
same scaling properties and regular behaviour. As an example, instead of eq. ( p.44| ) 
one could take 

1 1 

"^^^^ ^ 2iD(\-a) ^ |g|Z)(i-a) _^ 2£^(i-") ' (3.45) 

where I = 1/M. At small s (near the light cone) or large space/time scales, v — )■ 
const; at intermediate scales, v has the power-law behaviour which we will assume 
from now on. 

Without risk of confusion, we use the symbol d to denote total derivatives. 
Noting that 

d, = ^ds, □ = + ^^d, , (3.46) 

s s 



the inhomogeneous equation (|3.42|) reads 



{dl + ^^ds - m^) G{s) = 5i^\x) , (3.47) 

where we added a mass term, [m] = 1, > 0. We first consider the Euclidean 
propagator and denote with r = a/x^x* + the Wick- rotated Lorentz invariant. 
In the massless case, the solution of the homogeneous equation is 



Da 



G = Cir'y-— , (3.48) 
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where C is a normalization constant. The right-hand side of eq. (|3.47 ) is not the 6 
defined in radial coordinates. To find the latter, one notices that 



1 = j d''xv5^^\x) 
dr 5{r) , 

where = 2n^^'^ /T(D /2) is the volume of the unit D-ball. Therefore, 



dg{x)6i^\x) = / dg{x 



■5{r) 



(3.49) 



Hence, to find the propagator also for r = one can take some test function ip and 
compute 

/•+00 

VtoilCG,^) = \imVLD dr }CG(r)if(r) , (3.50) 

^-^0 Je 



where 



= dr (r^^-^dr) - r^^-'m' . (3.51) 



Therefore, 



/ + 00 
drG{r)dr [r^'^-^dr^) 

= GnD{2 - DaMO) , 

where we have used eq. ( |3.48| ) and integrated by parts once (boundary terms vanish). 
The last line must be equal to {S,ip), thus fixing G. Then, the Green function for 
m = reads 

This result enjoys several consistency checks. When a = 1, it is the usual Green 
function Go for the Laplacian in D dimensions with standard Lebesgue measure: 

r(f-i) 



limG(r) = Goir) 



In the UV limit a — )■ 2/D, one can expand eq. ( p.52| ) as r^/e = 1/e + Inr + 0(e). 
Up to a divergent constant, this is the logarithmic propagator in two dimensions, 
rescaled with a volume ratio: 

G,{r)= lim G{r) = ^G2{r) = ^hir. (3.53) 
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As a side remark, notice that G can be written as 



G{r) oc -^GD{r) , 

where 

l3 = D{l-a), 

and the fractional derivative can be defined as the Liouville derivative: 



dr-P^' 



r(7 + i) 
r(7 + ;5 + i) 



(3.54) 
(3.55) 

(3.56) 



The order of the derivative is negative, so that it is actually a fractional integral with 
clear meaning: Starting from the problem ^Gd = ^^^^ and inserting (heuristically) 
the identity d^d~^ , it replaces the second-order operator on Gd with a fractional 
differentiation of order 2 + /3 on G. We can conclude that 



G(r)ocG(i+^/2)(^) 



(3.57) 



i.e., the (massless) propagator is proportional to the Green function solving the 
pseudodifferential equation 



□ l+/3/2^(l+/3/2) ^ ^(D) 



which was calculated and discussed in 77, 78, 79, BO, Bl, 82] and reads 



1+ 



(3.58) 



(3.59) 



Indeed, after Wick rotation eqs. ( p.48| ) and ( p.59| ) agree up to the normalization. 
In this sense, in configuration space our field theory on a fractal is equivalent to a 
certain class of nonlocal models represented by eq. ( p.58| ) pT| , p2| , ^3] . 

We now consider the massive case (Helmholtz equation). The solution of the 
homogeneous equation CG = is 



G{r) 



(7) 



CiK Da_^{mr) + C9/0C. ^(mr) 



(3.60) 



where Ci^2 ai'e constants and K and I are the modified Bessel functions. Since for 
small m the solution must agree with the massless case ( |3.48| ), we can set G2 = 0. 
In fact, this is true only for a 7^ 2/D, as one can see from the asymptotic formulae 



KJz) 



Liz) 



-ln(f)-7 ifz/ = 



(3.61) 
(3.62) 
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where 7 is the Euler-Mascheroni constant. We shall discuss the case a = 2/D 
separately.^ 

To find the solution of the inhomogeneous equation, one exploits the fact that 
the mass term does not contribute near the origin. Expanding eq. ( |3.6(]| ) at mr ~ 
when Da>2 {C2 = 0), 



Da 
2 



which must coincide with eq. ( p.52| ) . This fixes the coefficient Ci and the propagator 
reads 

1 r(f) /mx^-i 



^W = -^ffAyl^j ^--i(-), (3.63) 

in agreement with the Helmholtz propagator (a = 1). 

Here an important remark is mandatory. In usual quantum field theory, the 
propagator G{r) is defined up to an immaterial constant C. By "immaterial" we 
obviously mean that □[G'(r) + C] = □G'(r). This is true in any dimension D and 
regardless the functional form of G{r), being it just a property of ordinary differen- 
tiation, nC = 0. Our model is still characterized by ordinary differential operators, 
so the set 

{G{r) + C \ C = const} (3.64) 

is an equivalence class defining the propagator. However, the nontrivial measure g 
will associate elements of this equivalence class with different Fourier transforms in 
momentum space. In the critical case a = 2/D, the difference will be in unusual 
terms _^ 

B{Ci, C2, C)-^ 

which resemble massless poles. If the theory is well defined, these terms will have 
to correspond to generalized Dirac distributions, and not to particle modes with 
arbitrarily chosen residue B{Ci, C2, C). Keeping this in mind will prevent us to fall 
into a false paradox when calculating the propagator in momentum space. 

Taking this issue on board, the case a = 2/D, mT^Ois straightforward: it is 
safe to just set a = 2/D in the noncritical propagator (|3.63|) , 

a(r) = --^Ko(mr). (3.65) 



®The comparison with field theory with fraetional powers of the d'Alembertian holds also when 
m 0. Recalling that 

one can continue this formula to any n and write the propagator as 
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Now to the Lorentzian theory. The fractal field models of |pT| and p2| both 



meet the Osterwalder-Schrader conditions. This encourages the expectation that a 
generic field theory on a Euclidean fractal, if well defined, should admit an analytic 
continuation to a theory in Lorentz spacetime. In fact, the Euclidean partition 
function eq. ( p.40D is a Schwinger function endowed with all the properties required by 
the Osterwalder-Schrader theorem: it is analytic, symmetric under the permutation 
of arguments, Euclidean covariant, and satisfies cluster decomposition and refiection 
positivity. Consequently, we can analytically continue the Helmholtz propagator 
( |3.63|) to the Klein-Gordon propagator according to the prescriptions: (i) multiply 
G times the imaginary unit i, due to Wick rotation of the time direction; (ii) replace 
with + ie, where the positive sign of the extra infinitesimal term corresponds 
to the causal Feynman propagator. Summarizing, 

z r(f) 



G{s) 



m 



Da 
4 



m 



while the massless propagator is 

G{s) 



s>0, (3.66) 



(3.67) 



nD{2 - Da) 

in accordance with eq. ( p. 591) . Equations (|3.53|) and ( p.65|) are continued similarly. 

The propagator for timelike intervals is just the analytic continuation of the 
former. In the massive case, it is proportional to the Hankel function of the first 
kind H^^^^^_^. One last thing to check is what happens on the light cone. The 
propagator for Da = 2, eq. ( |3.53| ), is (proportional to) the usual propagator in two 
dimensions, so nothing special occurs. Taking instead the definition ( |3.44 ), setting 
a = 1 in eq. (|3.67|) (contribution of m negligible), for even D = 4,6, . . . one obtains 
^ p. 94) 



G{s) 



oc 



2 r 



(D/2-2)! 



PV — - i7r6{s'^ 



(3.68) 



where PV denotes the principal value. In D = 4 this reduces to the Plemelji- 
Sokhotski formula. Translated into momentum space, the 6 states, as usual, that 
massless particles propagate at the speed of light on the light cone (Huygens' principle 
f|,|8|). 



3.3 Propagator: momentum space 

Since G has been argued to be proportional, in configuration space, to the Green 
function of another well-known problem (the functional inverse of a fractional power 
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of the d'Alembertian), we can already guess its pole structure in momentum space: 
in general, it will exhibit a branch cut with branch point at fc^ = — m^. 

To calculate the propagator in momentum space, we can start from the Euclidean 
one and then analytically continue the result as usual. In the Lorentzian propagators 
the substitution k"^ — )■ |kp — {k^Y — ie is understood. 

In the presence of a Lebesgue-Stieltjes measure, the Fourier transform must be 
modified so that it is consistent with the definition of the Dirac distribution eq. (|3.1|) . 
The Fourier-Stieltjes transform of a function G{x) and its inverse are defined as 
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(3,69) 

G{x) = [ dg{k) G{k) . (3.70) 



(27r) 



The measure in eq. ( p.70|) is such that momentum and configuration space have the 



same dimensionality. The Fourier-Stieltjes transform of a function G is the Fourier 
transform of vG: 

F^[G] = F[vG] = F[r^("-i)G'(r)] . (3.71) 
In particular, the Fourier-Stieltjes transform of 6^^^ is 1. When a 1 one has 



F41] = F[r^(°-i)] 

F (^) 1 



r 



D{l-a) 
2 



f.Da 



= i2nfSi''\k), (3.72) 

so Si^\ the source in eq. ( p.42| ), is a power-law distribution. 

Equation ( p.71| ) tells us the form of the massless propagator when a ^ 2/D 
(transform of eq. H: 

G(k) = D{a-1) 2^Da^ 

D-2 1 

(3.73) 



- 2 F ■ 
Notice that: 

• The (Lorentzian) propagator has a A;^ = pole in the (Re k^, Im k^) plane. From 
the world-fractal point of view (that is to say, looking at the pole structure of 
G{k) rather than of v{k)G{k)), the spectrum has the usual support at fc^ = 0. 

• When a = 1, one obtains G{k) = —l/k"^ and the free wave solution. 
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• For general a > 2/D, the sign of the residue is always negative, which ensures 
the absence of ghosts. However, its value is not 1 but given by a geometric 
factor. This is expected as the effective theory in the world-fractal is not unitary 
and some probability is exchanged with the D-dimensional topological bulk. 

In the critical case a = 2/D eq. ( |3.73D is ill-defined; this is not a problem, since one 
should start from eq. (|3.53|) : 

G^k) = J-Fp-^lnr] 

D \B-lnk'^ , , 

1 ^ (3-74) 



2 y fc2 

where B = iIj{D/2 — 1) — 7 + In 4 and ip is the digamma function. Some remarks: 

• For D 2, the Lorentzian propagator has a branch point at fc^ = 0. 

• In the limit D ^ 2, the propagator is G ~ —l/k"^. 

• The term B/k"^ does not represent a particle mode. In fact, it is nothing but 
the fractal Dirac distribution eq. ( ^.72| ) when a = 2/D. In other words, the 
Fourier-Stieltjes transform of the equivalence class ( p.64| ) is unique up to a 
Si^\k) term. The spectrum consists in a quasiparticle continuum of modes 
with momentum A;^ < 0. 

In eq. (|3.64|) there exists a particular element such that B = identically, which we 
typically call the propagator. Equation ( |3.74| ) is then replaced by the unambiguous 
expression 

D \lnfc2 
2 ~ ~P 



G*{k) = -{--l]^. (3.75) 



The massive case is slightly more complicated and as an exercise we will calculate it 
explicitly. The transform of the propagator in radial coordinates is 



G{k) = j dnDdrr^-^v{r)G{r)e-'''-\ 

The integrand is not radial but one can choose a frame where kf^x^ = —kr cos 6, 
k = I /c^ I , and the angular integral reads 

J dnDC-'^ '' = ^D-i de {sin 9)^-^6'''"'°''^ 

In the last line we used formula 3.915.5 of Then, 

G{k) = T (^^^ (^0 ' drr^'^-'^[nnG{r)]Jn_,{kr) . (3.76) 
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Now we take the massive propagator (|3.63|) (a 7^ 2//^): 

G,,)- r(f) ^2^«-V»^^-. 



/ drr~^ ^'^^Kdc ^(mr)jD i(kr) 
Jo ^ ''^ 



Da I - D. 

2 ' ' 2 ' m2 



(3.77) 



where we used formula 6.576.3 of ||85| and F is the hypergeometric function 



(3.78) 

1 ai \.{c + n)nl 

n=0 \ / \ / \ / 

When a = 1, one obtains the usual massive propagator in D dimensions (formula 
9.121.1 of [|5 



fc^ + m"^ 

In the limit m — (which does commute with the analytic continuation of eq. ( |3.77|) ), 
one recovers eq. ( |3.73| ) up to a term 0(m^°~^), which is negligible by virtue of the 
no-ghost condition a > 2/D. 

When a = 2/D, the transform of the critical propagator eq. ( p.65| ) is 

G*ik) = ^ 5 (3.80) 

§-iy-^ + 0{m'\nm'), (3.81) 

where in the second line we have dropped 0{k~'^) terms and considered the massless 
limit for comparison with eq. (|3.75|) . When D = 4, by virtue of formula 9.121.6 of 
BsH , eq. ( I3.80D gives exactly (i.e., no 5^ terms) 

\ ' 711? 



a(A:) = -^ln 1 + — ) , (3.82) 



which is precisely the generalization of eq. ( p.75p . G^{k) has a branch cut for Re k^ < 
— -\/|kp + w?. Therefore, the spectrum of the theory has a continuum of modes with 
rest mass > m. 

Interestingly, the convergence domain of eq. ( |3.77| ) gives a bound on the topo- 
logical dimension D. The propagator is a hypergeometric series with convergence in 
the unit circle \z\ < 1, where z = —k'^/m'^. Let 

Da D , , 

p^ — -- + l. (3.83) 

On the unit circle \z\ = 1, there are three cases: 
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• The series diverges ii p > 1. As a constraint on a it yields a > 1. This case is 
excluded by construction, although the limit case a = 1 is well-defined. 

• The series converges absolutely if p < 0, corresponding to a < 1 — 2/ D. 

• The series converges except atz = lifO<p<l, giving 1 — 2/D < a < 1. 

• When a = 2/D, the series diverges ii D < 2, converges absolutely if -D > 4, and 
converges except at z = 1 if 2 < D < A. These bounds are basically unchanged 
if one considers the analytic continuation of F to the massless case. 

The standard physical setting (convergence of the propagator in the unit disk and 
on its boundary except at the singularity on the real axis) is naturally recovered 
only for 



where we included both extrema of the interval by analytic continuation. 

Before concluding the section, we reconsider the issue of the superficial degree of 
divergence of Feynman graphs in the UV (see, e.g., p, ^ for an introduction to the 
subject and references). Consider a one-particle- irreducible subdiagram with L loops, 
I > L internal propagators and V vertices. The superficial degree of divergence 6 is 
the canonical dimension of all these contributions. 

Each loop integral gives [dg{k)] = Da, while the propagator, in any dimension 
and for any value of a, has [G] = —2. For the scalar field theory, interaction vertices 
do not carry dimensionality. Overall, 



When a = 1, one gets the standard result in D dimensions. In the critical case 
a = 2/D, 6 < and one has at most logarithmic divergences. When a < 2/D the 
theory is superrenormalizable. In the case of gravity also vertices contribute 0, |], 
each with a factor of 2 (number of derivatives). Then, 6 is bounded by the dimension 
of operators which already appear in the bare action. 

4. Gravity 

Having studied the properties of a scalar field on an effective fractal spacetime, we 
turn to gravity. Our conventions for the Levi-Civita connection, Riemann and Ricci 
tensors, and Ricci scalar are 



2< D<A, 



(3.84) 



S = L{Da - 2) - 2(/ - L) < L{Da - 2) . 



(3.85) 




(4.1) 
(4.2) 
(4.3) 
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The Ansatz for the gravitational action is 

Sa = ^2 j V^{R-2X- cod^vd^v) , (4.4) 

where g is the determinant of the dimensionless metric g^,^, = SttG is Newton's 
constant, A is a bare cosmological constant, and the term proportional to u has been 
added because v, like the other geometric field g^^, is now dynamical. The couplings 
have dimension 

[K'^] = 2-Da, [A] =2, [uj] = 2D{a-l) +A. (4.5) 

In spacetime with D = 2 topological dimensions and trivial measure weight v = 1, 
the Einstein-Hilbert action is a topological invariant and there are no dynamical 
degrees of freedom. This is not the case of eq. (|4.4| ). 

To describe the flow from the UV to the IR fixed point, we should add relevant 
operators also into the gravitational action. (The relevant operators in the matter 
sector are minimally coupled with gravity and they would not be enough.) This is 
done in the same way as for the matter sector, eqs. ( p.2ip and ( |2.22| ), with possibly 
different Mg ^ M^. The effective Newton constant then runs from a UV bare 
dimensionless constant to an IR value 

k',^ ~ 4vM^^ . (4.6) 

Note that Kjj^ is not necessarily the observed Newton constant k^^^, in which case 
Mg ~ mpi is the Planck mass. As one can see from the equations of motion, k^^^ will 
depend on the background as well as on the scale of the problem. 

Those in eq. (|2.21|) are not all possible relevant operators. Higher-order Riemann 
invariants Riem" ([Riem"] = 2n, n > 1) are irrelevant, but one could introduce 
lowest-order terms of the form (1 -|- Riem)" with n < 1. These terms might be 
expected in the context of fractal models, where fractional derivatives can (but not 
necessarily) find their natural setting. Considering all these operators, one ends up 
with a term 

^ Riem \ " 



dnc(n) 1 + , (4.7) 

\ mass^ / 

where cin) are arbitrary dimensionless coefficients. In the unrealistic case where 
they are all equal to 1 and Riem = i?, the integral can be summed explicitly to 
a nonpolynomial functional which admits Minkowski as a vacuum and yields the 
Einstein-Hilbert Lagrangian with cosmological constant at small R (in mass units): 

^ -l + ^-^ + 0(/?^). 



ln(l + i?) 2 12 



"'^^If they were regarded as necessary, all terms should make their appearance in any sector. 
This is a different definition of the theory which we will not consider. 
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This particular f{R) model might be of some cosmological interest. However, it is a 
toy model and we shall not continue its discussion. In fact, we are interested in the 
equations of motion near the UV fixed point, so we will ignore relevant operators 
from now on {Mg = = 0). 

4.1 Einstein equations 

Assuming that matter is minimally coupled with gravity, the total action is 

S = Sg + S^, (4.8) 

where Sg is eq. ( |4.4| ) and = J dg^^—gCrn is the matter action. The derivation of 
the Einstein equations is almost as in scalar-tensor models. We shall repeat it here 
to make the presentation self-contained. To find the equations of motion we need 
the variations 



SV^=-l9,uV^Sg^\ (4.9) 

6R = {R^, + g^uU- V,) 6g^^ , (4. 10) 

where VuV^ = dyV^ — TJ^^Va is the covariant derivative of a vector and the curved 
d'Alembertian on a scalar (p is 

□0 = ^9^(73^9^0) . (4.11) 



-9 

The Einstein equations SS/ 5g^^ = read 

= (4.12) 



^IJ,U = RfMU ~ -QtMuiR ~ 2 A) + 9fJ.U— ^ 



+UJ (^gf.udavd^v - d^vdyV^ , (4.13) 

T..^-^^--2'^^..X.,. (4.14) 

Taking the trace of eq. ( [4.12| ) gives 

-(^j-l^R + DX + {D- 1)^ cod^vd^v = (4.15) 
When taking into account the variation of the total action with respect to the scalar 

R-2X = -2k'^C^ - 00 {2vDv + d^vd'^v) , (4.16) 

eq. ( |4.15| ) becomes 

R+{D-l)^ + uj [DvDv + {D- l)d^vd''v] = (T/ - DC^) 



-2k' Ti^. (4.17) 
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Tu,y is the stress-energy tensor of matter. Its definition determines the continuity 



equation |8^. In fact, let 



SSr, 



(4.18) 



be the infinitesimal variation of the matter action with respect to the external fields 
6g^a^ and Sv. For the infinitesimal coordinate transformation ( p.l9| ), one has 



6v = a^dyV , 



(4.19) 
(4.20) 



where we used the definition of the Lie derivative for rank-2 and rank-0 tensors. 
Plugging eqs. ( [4.19| ) and ( [4.2(]| ) into ( [4.18| ) and integrating by parts, we get 



— / d^xa'^ 



(4.21) 



6Sra must vanish on shell (i.e., when the dynamical equations are satisfied). Using 
the properties of the Levi-Civita connection ( [4.1| ) and the definition of the covariant 
derivative of a rank-2 tensor, 

VT T^M _ f) n^fJ- _|_ T^o- _ -per rpfM 

V 9 



one finally obtains the continuity equation 



(4.22) 



which generalizes eq. (|3.20|) . 

If matter is a scalar field, it is straightforward to see that its equation of motion 
5S'm/50 = is eq. ( p. 41) with □ given by eq. ( [4.11|) , in agreement with eq. ( [4.22|) . 



The continuity and Einstein equations are not independent because of the con- 
tracted Bianchi identities 2V^R^u = ^ uR and eq. ([4.16|) . The divergence of {v times) 
eq. ( [4.12| ) correctly reproduces eq. ( |4.22| ). The check takes into account that in the 
absence of torsion the covariant derivative commutes on a scalar, [V^, V u\v = 0, 
while on a vector [V^,, Vy]V„ = R^^^'^Vr- 

4.2 Cosmology 

With the notable difference that matter is nonminimally coupled with the scalar f , 
the equations of motion are similar to those of Brans-Dicke theory [^, which is 
well constrained by large-scale observations |Q. This fact and the foreign physical 
setting lead to an altogether different dynamics. 
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In this section we specialize to a Friedmann-Robertson- Walker (FRW) line ele- 
ment 

ds^ = g^^dx^'dx" = -df + a{tfgijdx'dx^ , (4.23) 
where t is synchronous time, a{t) is the scale factor and 

dr*^ 

gijdx'dx^ = + r^rffi^ (4.24) 

is the line element of the maximally symmetric {D — l)-dimensional space S of 
constant sectional curvature K (equal to —1 for an open universe, for a flat universe 
and +1 for a closed universe with radius a). In four dimensions, dfl"^ = d9'^+sin^ Odip^. 
Quantities built up with the spatial metric gij will be decorated with a tilde. On this 
background, the only nonvanishing Levi-Civita and Ricci components are 



and 



V% = Hg,,, Vl^ = H5l ^{, = f:^, (4.25) 



Ro^ = -{D-l){H^ + H), (4.26) 

• ~ 2k 

R^i = Rij + [{D - l)H^ + H]gi,, R,j = —gij , (4.27) 

(X 

R={D-1)(^ + DH^ + 2Hj , (4.28) 



where 



H = - (4.29) 

a 

is the Hubble parameter (not to be confused with the Hamiltonian H of section |3.1| ) 
and we have exploited the symmetries of S . 

For simplicity we consider a perfect fluid (zero heat flow and anisotropic stress) 
as the only content of the universe: 

T^,u = {p + p) Uf^Uu +pgt,u, (4.30) 

where p = Tqq and p = */(/} — 1) are the energy density and pressure of the fluid 
and = (1,0, .. . , O)'^ is the unit timelike vector {u^u^ = —1) tangent to a fluid 
element's worldline. We also take a timelike fractal v = v{t). 
The 00 component of the Einstein equations ( [4.12 ) 



IS 



D A ^^i) 1 u .0 K A K 



— -1] H"^ + H- - -—^v' = — — p + — ^, (4.31) 

2 J V 2D-1 D-V D-l a?' ^ ' 

while combining that with the trace equation ( [4. 15] ) one obtains 

(4.32) 
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Other useful expressions can be found by suitable combinations of the dynamical 
equations. From eq. ([4.17|) , 

nv 

V 



R+{D-l)— + uj [DvUv -{D- l)v^] = -K^{p + p) , (4.33) 
and eq. ( ^4.16| ) one gets 



2A + - 1)— + - 2)co{vnv - v^) = -K^{p - p) , 



while from the 00 component of ([4.12|) and eq. ([4.16|) , 



H^ + H -H- + ^^vDv 
V D-1 



K 



ip + p) 



(4.34) 



(4.35) 



D-1 

If p + p 7^ 0, one can combine eqs. ( [4.33| ) and ( [4.35| ) to get a purely gravitational 
constraint: 



H+(D- 1)H^ +'^ + — + H-+ coivDv - v^) 

V V 



The continuity equation ( [4.22] ) contracted with —u'^ is 



P + 



{D-1)H + 



(P + P) = 0, 



0. (4.36) 



(4.37) 



and the Hubble 



where we used the definition of proper-time derivative, m^V^ ■ 
expansion = V^u^ = {D — 1)H in covariant formalism |8^, Note that this 

is not the volume expansion as in standard general relativity, as the latter is actually 

rfln(a^-if) 



e 



dt 



(4.38) 



which is the square bracket in eq. ( [4.37|) . For a barotropic fluid p = wp, the continuity 
equation is solved by p ~ {a^~^v)~^^'^'^\ up to some dimensionful prefactor. In 



general w is not a constant and we define the effective barotropic index 

^ (4.39) 

iV2 - V, 



The scalar field is a particular case of perfect fluid, with p = 



p = 0^/2 + V, and = -d^ 



. Equation ( [4.37| ) becomes eq. 



+ 



{D-l)H+- 

V 



+ 1/' = . 



(4.40) 



When V = 1 and = 4, we recover the standard Friedmann equations in four 
dimensions, eqs. ([4.31|) and ( |4.32|) (no gravitational constraint): 

^ A 

y 



K 



H' = -p + --- 

1 1 



3 a- 

.2 



X 



H' + H = --{3p + p) + -. 

3 



(4.41) 
(4.42) 
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On the other hand, for the measure weight 

v = t-^ , (4.43) 

where (3 is given by eq. ( p.55| ), the gravitational constraint is switched on. Then, the 
above equations should be taken cum grano salis. The UV regime, in fact, describes 
short scales at which inhomogeneities should play some role. If these are small, the 
modified Friedmann equations define a background for perturbations rather than a 
self-consistent dynamics. 

Modulo this caveat, we can look at flat (k = 0) background solutions in the deep 
UV regime with no cosmological constant. The gravitational constraint ( [4.36| ) is a 
Riccati equation in H which, together with the useful formulae 



V t ' V t 



(4.44) 



fixes almost completely the background expansion, regardless the matter content. A 
direct consequence of this overdetermination of the dynamics is that there are no 
vacuum solutions (p = = p). This might not happen for other fractal profiles than 
eq. ( [4.43| ), but at early times the measure must scale as eq. ( |4.43| ). This feature. 



therefore, is robust. 

Let us consider the cases w = and a; 7^ separately. When {3 = D — 2 = 2 
(UV regime), the a; = solution is 

a{t) = + , (4.45a) 

, _ H _ [f-{U + ^V2f2)c][f - (14 - 3v/22)c] 

(t9 - 2c)2 ' ^^■^^''^ 

where c is an integration constant. The energy density and pressure which solve all 
the equations simultaneously are 



3 [t9 + (v/IO-3)x/lOc][t9 + (v/IO + 3)v/IOc] 

These expressions are sufficient to characterize three cases: 

• c > 0: The scale factor decreases {H < 0) from t = until t = t^, = (2c) 
where the universe bounces (if* = 0, a* = 3^/^2~^/^c^/^). From t = to 
t = ti = [(14 + 3y22)c]l/^ the universe expands in superacceleration (e < 0), 
while for t > ti the expansion is only accelerated. The energy density p is 
negative for t > t^,, while the pressure p is always negative. 
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c = 0: Linear (decelerating) expansion, a = t, while p = p < always. 

c < 0: The universe expands in deceleration from a big bang event at t = 
to = |c|^/^ The energy density and pressure are negative for t > |4c|^/^ and 
t > |(v^ + 3)v^c|^/^ respectively. 



All these scenarios need a matter component with non-positive definite energy den- 
sity, so they are excluded if only ordinary matter is allowed. We envisage four simple 
modifications of this result. One is to change the flat prescription K = 0. Another is 
to consider the above formulae only asymptotically, since the simple measure profile 
eq. ([4.43|) is certainly valid only at early times. The only case where the universe 
expands at small t with p > is c < 0, for which one has superstiff matter {w{t) > 1). 
A third option is to allow for a nonzero geometric contribution U{v), a potential for 

V. 

A fourth possibility is that u ^ 0, which does lead to interesting cosmology. 
There is only one real solution to the gravitational constraint, namely, 

1/11 13 3uj^ 



2 220J HMlM (4 471 

where $ (also denoted as iFi or M) is Kummer's confluent hypergeometric function 
of the first kind: 

£j^£!. (,48) 

The expressions for p and p are 

_ 2(2u; + 3t4)(3u; + M^) $ f ; 



132tl0$fll. 13. 3a;N2 

24a;(2c^ + 3t4)$(^;f;|^) 



(4.49) 



2(u; + 3t^)(6a; + 5t^) 48.;^ $ f ; 

^ tlO ^ miO$fll. 13. 3a,N2 • ^^"^^^ 

V 4 ' 4 ' 2t4y 

We can use the asymptotic forms of $ to have some semi-analytic insight of the 
system. When z — )■ — cxd, 

, , ^^-oo r(6) ^^"^1(0 + ?^) r(6) 2;"" 

^> ~ f(irZ^(-^) ■■ + (-^) -g r(a) r(6-a-„) ^- <^'^i' 
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where is some finite order. On the other hand, 



ft; .) M e-.- + ,-,--^flit_l±A^K_ t^ . (4.52) 

r(a) ^ r(o — a) r(a — 77.) 

At late times the scale factor decreases and the fluid behaves effectively as phantom 
matter {w < —1): 

a ^ , Hr.-l, (4.53a) 



1 

'2' 



(4.53b) 



24 30 , , 

P ~ ^, P (4.53c) 

w ~ (4.53d) 

At early times we must distinguish between positive and negative u. For a; > 0, the 
universe is contracting and the fluid behaves like an effective cosmo logical constant: 

a ~ const X — , H ~ -— , (4.54a) 
e ~ (4.54b) 

P ~ ^ , P ^ , (4.54c) 

w 1. (4.54d) 

These quantities are plotted in figures |I] and || for = +1. 

For w < 0, at early times the universe expands and accelerates, even if the perfect 
fluid is stiff: 

a ~ const X , H (4.55a) 

06 



3 

5' 

2\uj\ 2\uj 



(4.55b) 



P ^ , P ~ , (4.55c) 

w ~ 1. (4.55d) 

These are plotted in figures ^ and |^ for uo = —1. 

The null energy condition is violated, so none of these scenarios can be realized 
by an ordinary scalar field, for which p + p = 0^ > 0. The system is invariant under 
time reversal, so the model with a; > 0, run backwards in time, also describes an 
expanding super accelerating (e < 0) universe filled with phantom matter. Contrary 
to standard general relativity, the energy density increases. This can be explained 
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1 2 3 401234 




012 3 401234 

t t 

Figure 1: The scale factor a, Hubble parameter slow-roll parameter e, energy density 
p (thick line) and pressm'e p (dashed line) for u = +1. 




12 3 4 

t 



Figure 2: The equation of state w = p/p ioi ui = +1. 

by recalling that the fractal geometry causes dissipation. The observer, in this case, 
experiences an incoming flux of energy from the four- dimensional bulk. 

The model with a; < is a universe which expands in acceleration. At some 
point the expansion is quasi de Sitter (notice the small plateau of H in figure |^), 
and after a brief period of superacceleration the cosmic expansion decelerates, stops, 
and reverts. Again, the fluid behaves unusually: during the accelerated expansion 
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Figure 3: The scale factor a, Hubble parameter H, slow-roll parameter e, energy density 
p (solid line) and pressure p (dashed line) for u = —1. 




0.5 1.0 1.5 2.0 2.5 3.0 

t 



Figure 4: The equation of state w = p/p for to = —1. The maximum, which lies beyond 
the frame, is w ^ 17.383 at t ^ 0.851. 



the equation of state is between dust and stiff matter (0 < w{t) < 1), while the 
energy density decreases. Then, after a non-monotonic transitory period around 
the inversion point, the equation of state violates first the strong and then the null 
energy condition. Eventually the energy density still decreases even during the cosmic 
contraction, because the world-fractal dissipates into the bulk. 
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The overall cosmological picture may be regarded as problematic for any value 
of uj. On one hand, there are no vacuum solutions. On the other hand, solutions 
with matter require fluids violating most or all energy conditions, thus signalling an 
unstable or unrealistic field (however, see below). The a; < case does so explicitly, 
since f is a ghost. 

This might be cured by allowing a nonzero intrinsic curvature, more complicated 
matter profiles, or a nontrivial potential for v}^ However, the most natural possibility 
is that a classical FRW background, either exact or linearly perturbed, is not realistic. 
Then, one would have to treat the UV limit as highly inhomogeneous. This is not 
at all unexpected, as we are dealing with quantum scales where the minisuperspace 
equations (maximal symmetry) are likely to fail. 



5. Conclusions and future developments 

There are several avenues of investigation left to explore. Here we mention just three. 

• Quantum field theory. Many aspects of the field theory have yet to be fully 
understood: among the most important are renormalization, the hierarchy 
problem and the physical significance of the UV propagator and the "natural" 
bound eq. ( p.84|) . Even in the Minkowski embedding, the fractal structure may 



have interesting properties we have not discussed here. For instance, because 
of violation of translation invariance the Fermi frame does not exist and, as in 
general relativity, parity is strictly a local symmetry. 

Other formulations of the theory including fractional derivatives would modify 
the propagator and dissipation properties and, pending a suitable definition 
of the kinetic terms, they could lead to a more transparent physics. Also, 
for simplicity we have defined an action on an embedding spacetime without 
boundaries, but a very interesting alternative is to consider scenarios with 
boundaries; for example, a field theory defined on MJ^ would be closer in spirit 
to the unilateral fractional and Weyl integrals of fractal classical mechanics, 
eq. (p^) with t — )■ +00. The propagator and several other features should 
change accordingly but, mutatis mutandis, the main idea of a fractal universe 
would still be valid and maintain the same motivations. 

Cosmology. If the system quickly flows to the IR fixed point, dissipative effects 
might be negligible on cosmological spacetime scales, with a notable exception. 
At late times, an imprint of the nontrivial short-scale geometry might survive 
as a cosmological constant of purely gravitational origin (pressureless matter 



^^Different geometric profiles would not work. For w > 0, the null energy condition is violated 



at early times {\t\ > 1 in the backwards model) where eq. (4.43) holds. For w < 0, w is a ghost by 
definition. 
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does not contribute to it). A detailed study may reveal whether the behaviour 
of the effective cosmological constant ( p.l5|) is compatible with observations. 

On the other hand, the UV regime may be relevant in the early universe, espe- 
cially during inflation. We have seen that the fiat background dynamics cannot 
be realized by a scalar field. It would be interesting to study the consequences 
of a nonvanishing intrinsic curvature or potential U{v). If, even in that case, 
ordinary matter (in particular, a scalar field) were not allowed, then one would 
have to abandon maximal symmetry and standard perturbative techniques of 
inflationary cosmology. However, an appealing alternative is to assume that 
matter is actually a condensate field stemming from a fermionic sector. It is 
known that a condensate violates the null energy condition, its mass-gap ef- 
fective energy density being negative in certain regimes 53 1. The physics 



of condensation is far from being exotic and is under good control. It would 
be interesting to see whether a Dirac sector with four-fermion interaction is 
renormalizable on a fractal and undergoes a condensation phase. 

Extra dimensions. In applications of the model we assumed that the topological 
dimension of embedding spacetime is D = 4. An interesting alternative is a 
universe with extra topological dimensions, D > 4. In this scenario the value 
of a would change and the IR limit should be realized in combination with a 
suitable compactification mechanism. 
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